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Abstract. In this presentation, I discussed a) the charm total cross-section and its 
comparisons to measurements at other beam energies and pQCD calculations; b) the 
scmileptonic decay of charmed hadrons and the sensitivity of non-photonic leptons to 
00 ■ charm quark collective flow and freeze-out; c) semileptonic decayed electron spectrum 

at high transverse momentum, its comparison to FONLL in p+p and d+Au collisions, 
and heavy-quark energy loss in Au+Au collisions. 
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In relativistic heavy-ion collisions, charm quarks are believed to be produced in the 
early stage via initial gluon fusion and their production cross-section can be evaluated 
using perturbative QCD pQ. Study of the N bin scaling properties of the charm total 
cross-section in p+p, d+Au and Au+Au collisions can test if heavy-flavor quarks, which 
are used as a probe, are produced exclusively at the initial impact. The interactions 
of heavy quarks with the medium provide a unique tool for probing the hot and dense 
matter created in ultra-relativistic heavy-ion collisions at the early times. At RHIC 



energies, heavy quark energy loss 0, charm quark coalescence jSl IH EI, the effect 
of J/ip production from charm quark coalescence on the interpretation of possible J/i/j 
suppression due to color screening [7j, and charm flow [HJ [HI| have been proposed 
as important tools in studying the properties of matter created in heavy ion collisions. 
The last three effects depend strongly on the charm total cross-section and spectrum at 
low p T . 

Since the beginning of RHIC operation, PHENIX and STAR collaborations have 
made pioneer measurements in charm related physics [TTJ H21 H3J 13 E3 EH] t- New 
measurements presented at this conference are: 

(i) muon spectra at forward rapidity (1.4 < \y\ < 2.2) from charm semileptonic decay 
by PHENIX Collaboration [IS] . This enables us to study the rapidity dependence 
of nuclear effects of charm production. 

(ii) muon spectra at low pt (0.17 < pr < 0.25 GeV/c) from charm semileptonic 
decay by STAR Collaboration 15J. This improves the charm total cross-section 
measurements and better constrains the charm spectrum for studying the charm 
radial flow. 

| The overviews of charm elliptic flow and quarkonium can be found elsewhere |17l 118) 
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Figure 1. Total cc cross-section per nucleon-nucleon collision vs. the collision energy 
(•V^iviv)- The dot-dashed and dashed lines are default NLO pQCD and PYTHIA 
calculations. The solid line is NLO pQCD calculations with /i^ changed from 2m c 
(dot-dashed line) to m c . 



1. Charm total cross-section 

It is difficult to directly reconstruct charmed hadrons and single electrons from charm 
semileptonic decay in hadron-hadron collisions with high precision at low p?, where 
the yield accounts for a large fraction (~ 85%) of the total cross section [TTJ [JJJJ 
12J. The difficulties are due to short decay distance (cr ~ 100/xm) and large 
combinatorial backgrounds in charmed hadron decay channels, and the overwhelming 
photon conversions in the detector material, and 7r° Dalitz decays in electron detection. 
Nevertheless, the charm total cross-sections have been measured in d+Au collisions at 
RHIC by a combination of the directly reconstructed low px D° — > Kit and the non- 
photonic electron spectra [XI], and by non-photonic electron spectra alone at pt > 0.8 
GeV/c jTSJ HHj. Although the systematic and statistical errors are large, the result 
indicates a much larger charm yield than predicted by pQCD calculations [TTJ UJ . It was 
argued that results with small renormalization (fragmentation) scales shown as solid line 
in Fig. [Hare not reliable calculations [JJ. A new method [2Hj was proposed to extract the 
charm total cross-section by measuring muons from charmed hadron semileptonic decay 
at low pt (e.g. 0.16^;>r^0.26 GeV/c). Since muons in this pt range are a very uniform 
sample of the whole charmed hadron spectrum, the inferred charm total cross-section is 
insensitive to the detail of the charm spectrum. 

1.1. Are we (RHIC) consistent? 

PHENIX and STAR collaborations have used several methods and techniques to extract 
total charm cross section in p+p, d+Au and Au+Au collisions [TTJ H2J EH H3J [JJ3 fTHj . 
The total cross-section measured at y/s NN = 200 GeV can be summarized as Fig. UJ and 
Fig. [21 In general, the measurements have smaller errors in Au+Au collisions than in 
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p+p or d+Au collisions. The agreements between PEHNIX and STAR are better in 
light systems than in central Au+Au collisions. The non-photonic electron spectra by 
PHENIX have smaller systematic errors than the corresponding STAR measurements 
while STAR Collaboration have two additional measurements from direct charmed 
hadron reconstruction and low-momentum muon spectra. These reflect the strengthes 
of the detectors, accordingly. There is substantial discrepancy of charm total cross- 
section between those extracted from PHENIX's non-photonic electron spectra and 
the combined (hadronic and semileptonic) fit results from STAR's measurements. The 
results in p+p and d+Au collisions show consistency within the errors. However, the 
discrepancy is about a factor of 2 in central Au+Au collisions while both have errors 
at 20% level. Part of the discrepancy may be explained by the different coverage of 
the two experiments [20]. Measurements from PHENIX non-photonic electron spectra 
cover < 15% of the dN/dy of the charm yields. Fig.Elright panel illustrates the possible 
difference between measurements at low pt and higher p?- It shows the dependence of 
the muon yield on power-law parameter (px) for a fixed total charm yield (details at 
Ref. [201 )• The yield is normalized to yields at (pt) = 1.3 GeV/c. Fig. [2] demonstrates 
that over a wide range in (pr), the muon yield is within ±15%. This is in contrast to the 
large variation of the electron yield integrated above px of 1.0 GeV/c, where a factor of 
8 variation is seen in Fig. [21 On the other hand, the large discrepance is difficult to be 
accounted for within reasonable range of parameters (e.g. (pt) need change from 1.4 to 
0.9 GeV/c for a factor of 2 change in extrapolated total yield) when both experimental 
errors are small in central Au+Au collisions. In addition, no obvious contradiction 
was observed from the electron spectra themselves and those are shown in Fig. 3 and 
Fig. 4 (details in later sections). With upgrades by PHENIX and STAR collaborations 
targeting at drastic improvement of precise secondary vertex detectors and continuous 
improvements of statistics and reduction/understanding of systematical errors, we may 
be able to better understand the difference in near future |2"T] . 

1.2. Are we alone? What's right and what's wrong? 

It has been shown that the charm cross-sections in hadron collisions at lower energy 
have large errors and in many cases are inconsistent. We examined the experimental 
data and how the total cross-section was extracted. We found that in most of the cases, 
the charm total cross-sections were inferred from measurements covering small phase 
space using formulae to extrapolate to the full phase space. In some cases, the functions 
used for extrapolation were not consistent and resulted in large discrepancy. In addition, 
the systematic error from the extrapolation has not been properly implemented. In X. 
Dong's thesis 1221, we tabulated all the measurements and got rid of those with large 
extrapolation or inferred from correlations. Although it eliminated many measurements, 
those kept show better consistency and were listed in ref. [221 E|- Now we can compare 
the experimental data from low energy to high energy to the pQCD calculations. At 
RHIC energies, the data points are in general above the default pQCD predictions. The 
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Figure 2. Left panel: Differential cc cross-section per nucleon-nucleon collision 
(dN/dy) as function of the number of collision participant nucleons (N part ). The 
solid line is default NLO pQCD and FNOLL calculations. The shaded band indicates 
uncertainty of the predictions. The data points are from PHENIX and STAR 
collaborations. Right panel: Lepton yields relative to the fixed total charm cross- 
section as function of power-law parameters (pr) for a charmed hadron transverse 
momentum spectrum. Solid line shows muon yields with a kinematics selection 
0.16 < pr < 0.26 GeVc and \yi\ < 0.5. Dashed line shows electron yields with pr > 1.0 
GeV/c. 

usual wisdom is that there is a factor of 2 uncertainty from pQCD calculations by varying 
the renormalization and fragmentation scales within a factor of 2. We noted that cross 
sections inferred from cosmic ray measurements are also much higher than the pQCD 
predictions. The cosmic ray showers required large charm total cross-section to account 
for the large muon yields and electromagnetic showers within the cosmic ray shower 
with incident nucleon energy at a few teens of TeV. On the other hand, the results at 
lower energy are consistent with pQCD calculations as shown in Fig. ^ and charmed 
hadrons at px > 5 GeV/c at Tevtron is within a factor of 2 above pQCD calculation 
We have also compared the non-photonic electron spectra from ISR energies to RHIC 
and to Tevtron energies. We found that there is no obvious inconsistency among the 
electron spectra. The direct D° measurement for px > 5 GeV/c in p + p collisions by 
CDF Collaboration was fitted to a power-law function and the semileptonic decayed 
electron spectrum was from the extrapolated spectrum. Detail of this study has been 
shown in Ref. |22H23|. 

2. Are leptonic spectra sensitive to charm flow? 

We propose a new method to extract the charm total cross-section by measuring muons 
from charmed hadron semileptonic decay at low px (e.g. 0.16<pt^0.26 GeV/c). Since 
muons in this px range are a very uniform sample of the whole charmed hadron spectrum, 
the inferred charm total cross-section is insensitive to the detail of the charm spectrum. 
Once the cross-section is determined, the electron spectrum at higher px can be used 
to sensitively infer the charmed hadron spectral shape. Meanwhile, we survey the form 
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factors used in charm semileptonic decays generated from Particle Data Group [21], in 
the PYTHIA event generator [TT) IT2"1 IT3*1 123] . by pQCD predictions pQ and from the 
CLEO inclusive measurement 26J. We find that the lepton spectra from these different 
form factors can be different by a factor of 1.5. 

2. 1 . Form factors 

The spectrum generated by the PDG is according to the form factor of charmed meson 
decays to pseudoscalar K+l+i/, vector meson K* + 1 + 1/ and non- resonance (Kii) + l + v 
where the K* mass is used for the (Kir) system. Since PYTHIA uses a simplified vector 
meson decay form factor [22], it tends to produce a softer electron spectrum. Both the 
parameterization by Cacciari pQ and formulae from the PDG agree with CLEO's electron 
spectrum j2SJ- In addition, we also find that although the charmed mesons (D ± and 
D°) from ^(3770) decay have a momentum of 244 MeV/c only and without correction 
of final state radiation [2E], it affects slightly its subsequent electron spectrum. 

2.2. Freeze- out 

The semileptonic decay greatly smears the spectrum and reduces the difference between 
the different spectrum shapes. However, it is clear that a reasonably realistic blast-wave 
parameterization of charmed mesons in Au+Au collisions is very different from that in 
d+Au collisions. There is also a significant difference between spectra with different flow 
(blast-wave function) parameters at 0.5 < < 1.5 GeV/c. Fig. 6 in Ref. [20J shows that 
there is a factor of 3 difference at pr = 1.5 GeV/c between late freeze-out (Tf a = 100 
MeV, f3 max = 0.9) 27 and early freeze-out (T fo = 160 MeV, (3 max = 0.6) [2E]. Current 
measurements of non-photonic electron spectra and direct charm spectra seem to be 
consistent with a decreasing trend even at pr — 1.0 GeV/c [H3III], which is likely due 
to multiple collisions and thermalization at low px with early freeze-out jUj and not due 
to pQCD energy loss. However, since the overall normalization is not known, an early 
freeze-out scenario can be interpreted as suppression of the charm total cross-section as 
well. This ambiguity can be resolved by a measurement of the total cross-section |20j. 
In fact, with the measurements by PHENIX and STAR collaborations presented in 
this conference, it was demonstrated that we are able to obtain freeze-out parameters 
based on a blast-wave assumption with reasonable errors [15j (also shown as solid line in 
Fig. 4). However, the errors on the current measurements are large in this pt range. We 
advocate improving the measurements of electrons in this pt range to assess if charm 
thermalizes in the medium jH]. 

3. Does it make senses: Color and Flavor? 

Recently, the measurements of high p? electrons from heavy-flavor semileptonic 
decays have posed challenges to our understanding of partonic energy loss in the 
medium [TJ1 03 El EE 122] ■ 
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Figure 3. Ratio of non-photonic electrons to semileptonic decayed electrons from 
charmed and bottomed hadrons in FONLL as function of pt- The data points are 
from PHENIX and STAR collaborations. 

3.1. p+p /FONLL is 5.5 at highp T ? 

It has been perceived that pQCD can reasonably calculate charm production because 
its quark mass is much larger than the non-perturbative scale. However, we have shown 
that the total cross section seems to be above pQCD calculations in several cases. On the 
other hand, the agreement should be better at higher pj>. Fig.OHsliows the measurements 
of non-photonic electron spectra and muon spectra in p+p collisions divided by the First- 
Order-Next-to-Leading-Log (FONLL) calculations jT] as function of p?- The pQCD 
calculations have a factor of 2 uncertainty. Although the experimental data have large 
errors at low px, it is in general larger than the prediction. In particular, the STAR 
data from Electromagnetic calorimeter and the PHENIX data from muon detector are a 
factor of 5 above the FONLL calculations. This presents a puzzling issue since charged 
hadron and pions [3UJ EI] at high p T can be quite well reproduced by NLO pQCD 
calculations. 

3.2. Color and flavor dependence of energy loss 

Let us ignore the discrepancy between data and pQCD in p+p collisions and study 
the nuclear modification function (Raa) i n Au+Au collisions. PHENIX and STAR 
collaborations has shown that Raa of non-photonic electron (presumedly from charm 
and bottom decays) is much smaller than 1 for pt^2 GeV/c shown in Fig. EJ In fact, 
it is very similar to that of pions with Raa — 0.2. Jet-quenching models incorporating 
collisional and radiation energy loss can not account for such large suppression [29J. 
This may be due to the large contribution of bottomed hadrons to the non-photonic 
electron in the model calculations, which assume that charm and bottom production 
scale the same way from pQCD to match data in p+p collisions [32] . 

If we look at the Raa at quark level for gluons, light quarks, and charm quarks, it is 



Overview of charm production at RHIC 



7 



1 1 1 1 — 

PHENIX prel.(0-10%) 
STAR EMC prei.(0-5%) 
STAR TOF prel. (0-1 2%/dAu) 
c-> e Late Freeze-out 
c-> e Early Freeze-out (large flow) 
c— > e Early Freeze-out (small flow) 




Transverse momentum p T (GeV/c) 



Figure 4. The nuclear modifiction function (Raa)o{ non-photonic electron spectra vs. 
transverse momentum. Curves show different freeze-out assumptions in a blast-wave 
model. The data points are from PHENIX and STAR collaborations. 



obvious that the difference between light quarks and charm quarks is quite small (much 
less than a factor of two at low p? and similar at high p?) calculated by Ref . [21] . The 
largest difference is between light quarks and gluons. This means that the deadcone 
effect of heavy quarks traversing the medium has much smaller effect on the nuclear 
modification factor than the color-charge factor of 9/4 between gluons and quarks. To 
study this, it is important to find differential experimental probes sensitive to gluon 
and light quark energy loss. It has been proposed [33J that anti-protons at high p T 
(10 GeV/c) in A+A collisions are mainly from gluon fragmentation while pions are 
dominated by light quark fragmentation. By measuring Raa-, P/p and p/ir ratios, we 
will be able to access the difference between gluon and quark energy loss. This was 
done by STAR Collaboration in a recent publication [HI] using the relativistic rise of 
ionization energy loss in TPC to separate protons and pions [33]. The results show 
that there is no difference between proton and pion Raa (or Rcp) and little centrality 
dependence of p/n ratio has been found [3T] . 

We now witness a set of data showing that the nuclear modification factors are 
the same among hadrons that are presumedly fragments from separated gluons, light 
quarks and heavy quarks. This apparently contradicts the prediction from jet quenching 
models, which successfully explain the light hadron production at high pt and dijet 
correlations [3U]- R doesn't necessarily mean that the general framework of jet quenching 
is invalid. This may imply that interesting phenomena are present besides the general 
idea of energy loss of energetic partons traversing dense medium. 

Experimentally, we need to measure the charm total cross-section with energy scan 
to map out its energy dependence and help constrain the pQCD calculations and charm 
coalescence into J/^fr. We also will be able to study radial and elliptic flows of heavy 
flavors with upgrades by PHENIX and STAR collaborations. We need to separately 
measure the nuclear modification factors of charmed and bottomed hadrons and possible 
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heavy-quark tagged jets to study the charm and bottom quark energy loss. 
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